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Abstract 



Recently it has been suggested to use the perigee of the proposed 
LARES/ WEBER-SAT sateUite in order to measure the secular pre- 
cession which would be induced on such Keplerian orbital element by 
a weak-field modification of gravity occurring in some Brane- World 
scenarios put forth by Dvali, Gabadadze and Porrati. This precession, 
derived for the first time by Lue and Starkman, amounts to ~ =f4x f 0~'^ 
milliarcseconds per year. In this paper we show that, according to the 
recently released EIGEN-CGOIC Earth gravity model which combines 
data from the dedicated CHAMP and GRACE missions, the quite 
larger systematic errors due to the Newtonian part of the terrestrial 
gravitational potential would vanish any attempts to detect a so small 
effect. Improvements in our knowledge of the Earth's gravitational 
field of, e.g., up to four orders of magnitude in the low degree even 
zonal harmonics would be required. Also the obtainable observational 
accuracy in reconstructing the LARES orbit rules out the possibility 
of measuring the LS effect with the perigee of LARES. The situa- 
tion is much more favorable in the Solar System scenario. Indeed, 
the non-Newtonian perihelion advance of Mars, which is the currently 
best tracked planet, has recently been measured with an accuracy of 
1 X 10^^ milliarcseconds over a 92-years time span; the LS precession 
is ~ =f4 X 10^^ milliarcseconds per century. A suitable combination 
of the perihelia and the nodes of Mars and Mercury, which disentan- 
gles the LS effect from the competing larger Newtonian and general 
relativistic precessions, might allow to reach a I-sigma 41% level of 
accuracy. 

Keywords: LARES, even zonal harmonics, planets of the Solar System, 
orbital motion, Brane- World scenario 
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1 Introduction 



The recently observed acceleration of the Universe could be explained either 
by invoking a small vacuum energy (or some effective vacuum energy) or 
by assuming that the Einstein gravity undergoes long-range modifications 
imitating a small cosmological constant. The second approach seems to be 
more promising because it is perturbatively stable under quantum correc- 
tions. In it the conventional gravity breaks down at very large distances, 
beyond a certain threshold Tc which is perturbatively stable and is fixed by 
the observations to a value ~ 10^^ m. In such a framework there is a class 
of Brane- World theories which, for r « r^, lead to Einstein-Newton grav- 
ity up to small corrections jjj. As a consequence of such tiny deviations, 
a secular advance of the argument of pericentre w of a test particle freely 
orbiting a central mass M occurs. 

According to the Dvali-Gabadadze-Porrati (DVP) model 01, which is a 
Brane- World theory with an infinite- volume flat extra spatial dimension, the 
analytical expression of the pericentre advance over one orbit of one object 
moving at distance r from a central body of mass M is j3] 

37r\/2 r3/2 , , . _i 

Aw = =F rad orbit , (1) 

4 

where Tc = 6 Gpci= 1.8514068 x lO^^ km and rg = 2GM/c^. The Lue- 
Starkman (LS) secular rate for any object orbiting a central mass along a 
circular orbit amounts to ~ =f4 x 10^^ milliarcseconds per year (mas yr~^), 
or ^4 X 10~^ arcseconds per century (" cy~^), only^. 

The possibility of testing such an effect in the Solar System scenario has 
preliminarily been investigated for the first time in |^. In a recent paper |S] 
it has been suggested, among other things, to use the perigee of the proposed 
Earth artificial satellite LARES in order to detect the LS effect. 

In the present work we analyze in more details such possibilities yielding 
a realistic evaluation of the error budgets by also including many relevant 
sources of systematic errors. 

4 parsec (pc)=3.085678 x 10^^ km. The value = 6 Gpc is quoted in 
^Note that in ~ 5 Gpc is assumed from fitting of both the Cosmic Microwave 
Background (CMB) and the type lA supernovas (Starkman, G. 2005; personal communi- 
cation): this yields 5 x 10~* " cy~^ for the magnitude of the LS effect. Its sign is affected 
by the cosmological branch 0: it is negative for the standard cosmological branch, and 
positive for the self-accelerated one. 
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Table 1: Orbital parameters of LAGEOS, LAGEOS II and LARES. 



Orbital element 


LAGEOS 


LAGEOS II 


LARES 


semimajor axis a (km) 


1.2270 X 10'' 


1.2163 X 10^ 


1.2270 X 10^ 


eccentricity e 


4.5 X IQ-^ 


1.4 X 10-2 


4.0 X 10^2 


inclination i (deg) 


1.10 X 102 


5.265 X 10^ 


7.0 X 10^ 



2 The Earth-LARES system 

The Earth artificial satellite L ARES/WEBER-S AT 0111 was proposed in 
order to measure the post-Newtonian general relativistic Lense-Thirring ef- 
fect on the orbit of a test particle to a high accuracy level (~ 1%) in the 
gravitational field of the Earth by suitably combining its data with those of 
the existing geodetic SLR (Satellite Laser Ranging) satellites LAGEOS and 
LAGEOS II. Their orbital parameters are in Tabled 

2.1 The observational sensitivity 

From Ar = eaAco 8 , it can be obtained that the accuracy in measuring the 
perigee precession over a given observational time span can be expressed 
as 5uj = brjea. By assuming a root-mean-square (rms) error of 1 mm 
in reconstructing the LARES orbit over, say, one year for a given set of 
dynamical force models one gets Jo; = 4 x 10"^ mas. 

With^ a = 3.6 X 10^ km and e = 2.8 x 10^^, as also suggested in [S], and 
by assuming a rms error of 1 cm, the accuracy in the perigee would amount 
to 2 X 10"^ mas. Note that the adopted values for the obtainable accuracies 
in r are optimistic; for example, the mm accuracy has not yet been fully 
achieved for the existing LAGEOS satellites. 

2.2 The systematic errors of gravitational origin 

The perigee of an Earth artificial satellite is affected by various kinds of 
long-period (i.e. averaged over one orbital revolution) orbital perturbations 
induced by the Newtonian Earth's geopotential ^01 • The most insidious 
ones are the secular rates induced by the even {I = 2,4,6...) zonal (m = 0) 
harmonic coefficients Ji of the geopotential which account for the departure 
of the Earth from an exact spherical shape. Their explicit expressions up to 

^It may be interesting to note that, with such an orbital configuration and i = 6.34x lO'^ 
deg the gravitomagnetic Lense-Thirring precession [?] would amount to -1.8 mas yr~^ only. 



4 



degree £ = 20 can be found, e.g., in The largest precession is induced 
by the Earth's quadrupole mass moment J2. For a moderate eccentricity its 
analytical expression is 

cLu 3 [RV _hf^ 5 . 2 .\ 

— = -ra — ^-77 2 sm I \ , (2) 

where n = \J GM /a^ is the unperturbed Keplerian mean motion and R is 
the Earth's equatorial mean radius. 

In |S] it is proposed to launch LARES in the so called frozen-perigee 
orbit characterized by the critical value of the inclination, i = 6.34 x 10^ 
deg, for which the J2-precession of eq.® vanishes. Moreover, it seems that 
Ciufolini suggests to use the LARES data together with those from LAGEOS 
and LAGEOS II in order to measure the Lense-Thirring as well. 

We will now show that such proposals are unfeasible. 

2.2.1 The impact of the even zonal harmonics of the geopotential 

Indeed, apart from the fact that the unavoidable orbital injection errors 
would prevent to exactly insert LARES in orbit with the required inclina- 
tion, it turns out that the impact of the other uncancelled precessions in- 
duced by the even zonal harmonics of higher degree, along with their secular 
variations, would totally swamp the LS effect. In Table |21 we use the cali- 
brated sigmas of the even zonal harmonics of the recently released combined 
CHAMP-hGRACE-Fterrestrial gravity EIGEN-CGOIC Earth gravity model 
[121 in order to calculate the mismodelled residual classical precessions on 
the perigee of LARES by assuming i = 6.34 x 10^ deg. It can be easily 
seen that the mismodelled precessions induced by the first seven even zonal 
harmonics are larger than the LS precession. Due to the extreme smallness 
of such an effect, it is really unlikely that the forthcoming Earth gravity 
models from CHAMP and GRACE will substantially change the situation. 
This rules out the possibility of using only the perigee of LARES. 

2.2.2 The linear combination approach 

In regard to the possibility of suitably combining the Keplerian orbital 
elements of the existing LAGEOS satellites and of the proposed LARES 
|14| I15| I16j in order to reduce the impact of the systematic errors of the 
geopotential on the proposed measurement, it is unfeasible as well. 

Indeed, the perigee is also affected, among other things, by the post- 
Newtonian general relativistic Einstein precession , whose nominal value 
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Table 2: Mismodelled classical secular precessions, in mas yr~ , of the 
perigee of an Earth satellite with semimajor axis a = 1.2270 x 10^ km, 
inclination i = 6.34 x 10^ deg, eccentricity e = 4.0 x 10~^ according to 
the combined CHAMP+GRACE+terrestrial gravity EIGEN-CGOIC Earth 
gravity model up to degree £ = 20 (1-sigma). For the precessions induced 
by the secular variations of the even zonal harmonics, referred to an obser- 
vational time span of one year, the values for J^, £ = 2, 4, 6 of ^H] have been 
used. Such precessions grow linearly in time. In the last three rows a - has 
been inserted because the corresponding mismodelled precessions amount to 
10~^ mas yr~^. The LS secular precession amounts to =f4 x 10~^ mas yr"'^ 
for an object orbiting a central mass along a circular path. 



Even zonal harmonic 


Mismodelled precessions (mas yr ^) 


J2 


1.26 X 10^^ 


j2 


4 X lO^'^ mas yr^^ 


■h 


4.390 


Ja 


1.406 mas yr^^ 




1.355 




6.63 X 10~^ mas yr~^ 


Js 


2.08 X 10^1 


Jio 


2.5 X 10-2 


Jl2 


2.5 X 10-2 


Jl4 


7 X 10-=^ 


Jw 




Jl8 




J20 
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for LARES is 3.280136 x 10^ mas yr~^, and by the Lense-Thirring effect, 
wliich, for i = 6.34 x 10^ deg, nominally amounts to -4.1466 x 10^ mas yr~^. 
This means that if we want to measure the LS precession independently of 
such quite larger Newtonian and post-Newtonian effects we would need ten 
Keplerian orbital elements in order to write down a linear system of ten 
equations in ten unknowns (the first seven even zonal harmonics, the LS 
effect and the two relativistic precessions) and solve it for the LS precession. 
Instead, we would have at our disposal, in principle, only four reliable Keple- 
rian orbital elements^: the longitudes of the ascending node of LAGEOS, 
LAGEOS II and LARES-which are affected by the Earth's geopotential and 
by the Lense-Thirring effect but not by the LS force-and the perigee of 
LARES. They would only allow to cancel out the general relativistic effects 
and just one even zonal harmonic. 

The data from an hypothetical satellite with a = 3.6 x 10"^ km, as pro- 
posed in [3], could not be used for the following reasons. If, on the one hand, 
the classical geopotential precessions would be smaller than the LS effect, 
apart from those induced by J4 and Jq according to EIGEN-CGOIC, on the 
other hand, the gravitomagnetic Lense-Thirring and the gravitoelectric Ein- 
stein post-Newtonian precessions on the perigee would amount to -1.8512 
mas yr (for i = 6.34 x 10^ deg) and 2.409956 x 10^ mas yr , respectively. 
This means that, even for such a higher altitude, one could not analyze only 
the perigee of LARES whose Keplerian orbital elements should be, instead, 
combined with those of the existing LAGEOS and LAGEOS II: with the 
nodes of the LARES, LAGEOS and LAGEOS II and the perigee of LARES 
it would be possible, in principle, to disentangle the LS effect from the post- 
Newtonian precessions and one classical even zonal harmonic. For such a 
high altitude the period of the node of LARES would amount to ~ 10^ 
days, i.e. tens of years. The tesseral Ki tidal perturbation, which is one 
of the most powerful harmonic time-dependent perturbations which are not 
cancelled out by the linear combination approach, has just the period of the 
node. Then, it would act as a superimposed linear bias over an observational 
time span of a few years. Moreover, as shown in jl81ll9j . when high-altitude 
satellites are included in linear combinations involving also lower satellites 
as the existing LAGEOS and LAGEOS II it turns out that the orbital ele- 

*Note that the other routinely and accurately laser-tracked SLR satellites which could, 
in principle, be considered are Ajsai, Starlette and Stella: the useful orbital elements are 
their nodes and the perigee of Starlette. However, since they orbit at much lower altitudes 
than the LAGEOS satellites they would practically be useless. Indeed, they are sensitive 
to much more even zonal harmonics of the geopotential so that they would greatly increase 
the systematic error induced by them. 
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merits of the higher SLR targets enter the combination with huge coefficients 
which amphfy all the uncancelled orbital perturbations. This would also be 
the case for the Ki tide affecting LARES. 



2.2.3 The impact of the odd zonal harmonics 

The perigee of an Earth artificial satellite is also affected by long-period 
harmonic perturbations induced by the odd {i = 3,5,7...) zonal (m = 0) 
harmonics of the geopotential. The largest perturbation is induced by J3: 
it has a sinusoidal signature with the period of the perigee. Its analytic 
expression is 



+ 



3 / R\ J3 sin a; 
2 V / esini(l — e^)^ 



35 2 2 
1 sin i cos i 



- sm I 
4 



1 sin i+ 



(3) 



For i = 6.34 x 10^ deg the period of the perigee of LARES, given by eq.©, 
is of the order of 10^ days; moreover, the second term of the right-hand-side 
of eq.© does not vanish. This means that the perigee of LARES in the 
critical inclination would be affected by an additional semisecular bias due 
to J3 which, over an observational time span of some years, would resemble 
a superimposed linear trend. According to EIGEN-CGOIC, its mismodelled 
effect would be < 5 x 10~^ mas yr~^. 

This additional bias should be accounted for both in the perigee-only 
scenario and in the linear combinations scenario. Note also that for a = 
3.6 X 10^ km and i = 6.34 x 10^ deg the period of the perigee would amount 
to ~ 10^ days, i.e. ~ 10^ years. 



3 The Earth-Moon system 

In ^ the possibility of analyzing the perigee of the Moon (a = 3.84 x 10^ km, 
e = 0.05490) is considered in view of the possible future improvements in the 
range residuals accuracy from the cm to the mm level in the the Lunar Laser 
Ranging (LLR) technique [25 • By assuming (5r ~ 1 mm, the uncertainty in 
determining the Moon perigee would amount to Jo; = 9 x 10~^ mas which is 
the magnitude of the LS shift over almost 2 years. The currently available 
accuracy in LLR range residuals is ~ 2 cm which maps into a 1.8 x 10^^ 
mas accuracy on the perigee. 
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Table 3: Orbital parameters of Mercury, Venus and Mars 
(http://nssdc.gsfc.nasa.gov/planetary/factsheet/). For the Astronomi- 
cal Unit (A.U.) we use the value 1 A.U.=149597870691 m. The angle e 
refers to the inclination to the ecliptic: the inclination i to the reference 
plane is almost e/2. 



Planet 


a (A.U.) 


e (deg) 


e 


Mercury 


3.8709893 x 10^^ 


7.00487 


2.0563069 x 10"^ 


Venus 


7.2333199 x lO^^ 


3.39471 


6.77323 X IQ-^ 


Mars 


1.52366231 


1.85061 


9.341233 X IQ-^ 



4 The Solar System scenario 

If we look at the planetary motions in the gravitational field of the Sun 
the situation is more favorable. The LS perihelion secular rates for all the 
planets, by neglecting the eccentricity of their orbits, amounts to^ =f4 x 10~^ 
" cy~^. It is undoubtedly a small value if compared, e.g, to the well known 
Einstein perihelion advance which, for Mercury, amounts to 4.29811 x 10^ " 
cy-^ 

4.1 The observational sensitivity 

Contrary to the Earth-LARES scenario, the braneworld effect lies at the 
edge of the present-day available orbit determination accuracy. 

The best tracked planet is Mars thanks to the various spacecraft mis- 
sions to it: Viking (1976-1982), Pathfinder landers (1997), Mars Global 
Surveyor (MGS) (1998-2003) and Odyssey (2002-2003). According to [12], 
the rms residuals of ranging for Viking are 8 m, for Pathfinder are 4.4 m, 
for MGS and Odyssey are 1.4 m. From the results of Table 4 of P2] for 
the nonsingular elements h and k, based on the processing of about 300000 
position observations (1911-2003) of different types used for the EPM2004 
ephemerides, it turns out that the formal standard deviation in the Martian 
perihelion advance is 1.5 x 10^^ mas; the braneworld shift over a 92-years 
time span amounts to 3.68 x 10^^ mas. 

The present-day situation for Mercury and, especially, Venus is less fa- 
vorable mainly due to the systematic errors induced by the mismodelling 

^Departures from such value due to the eccentricities are very hkely too smaU to be 
detectable, as can be inferred from the subsequent discussion. 
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in the planetary topography and the small value of the eccentricity of the 
venusian orbit. Indeed, from \22. the rms residuals of ranging from all ra- 
diometric data (1911-2003) amount to 1.4 km and 0.7 km for Mercury and 
Venus, respectively. The formal standard deviations in the perihelia advance 
are 7.6 x 10^^ mas and 8.772 mas, respectively. However, the future Messen- 
ger (launched in August 2004 by NASA) BepiColombo and Venus Express 
(ESA) missions should improve, among other things , also the accuracy of 
the ephemerides of Mercury and Venus. 

4.2 The systematic errors 

However, as in the case of the Earth-LARES scenario, the systematic errors 
of gravitational origin represent a very limiting factor. In Table 0] the nom- 
inal values of the major Newtonian and post-Newtonian general relativistic 
competing precessions are listed. It turns out that the accuracy with which 
the various competing effects^, are presently known would make difficult 
to use only one perihelion, apart from, perhaps, Mars. The major limiting 
effects are the classical N-body precession, the post-Newtonian general rela- 
tivistic gravitoelectric perihelion advance and the classical secular precession 
due to the Sun's even zonal harmonic J2. 

A very limiting systematic bias might be, in principle, that induced by 
the huge Newtonian N-body secular precessions whose nominal values are 
up to seven orders of magnitude larger than the effect of interest. In re- 
gard to their mismodelling, the major source of uncertainty is represented 
by the GM of the perturbing bodies among which Jupiter plays a domi- 
nant role, especially for Mars. According to j25j . the Jovian GM is known 
with a relative accuracy of 10~^; this would imply for the red planet a 
mismodelled precession induced by Jupiter of the order of 10~^ " cy~^. 
The GM of Saturn is known with a 1 x 10^^ relative accuracy |26j . How- 
ever, the ratio of the precession induced on the Mars perihelion by Sat- 
urn to that induced on the Mars perihelion by Jupiter is proportional to 
(Msat/-^jup)(ajup/«Sat)^ ~ 5 X 10~^. This would assure that also the effect 
of Saturn is of the order of 10^^ " cy~^. The situation with the precessions 
induced by Uranus and Neptune on the martian perihelion is even more 
favorable. Indeed, for Uranus (Mura/Afjup)(ajup/aura)^ ~ 1 x 10~^ and 
the relative uncertainty in the uranian GM is 2 x 10~^ -^o^ Neptune 
(MNep/Mjup)(ajup/aNop)'"3 x 10"^ and 6{GM)/GM = 2 x 10"^ In 
regard to the precessions induced by the most important asteroids, their 

®Note that the nominal values of the Lense-Thirring precessions are one order of magni- 
tude smaller than the LS effect, so that the gravitomagnetic force does not pose problems. 
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Table 4: Nominal values, in " cy~^, of the Newtonian and post-Newtonian 
general relativistic secular precessions affecting the perihelia and the nodes 
of the inner planets of the Solar System. For a given planet, the pre- 
cession labelled with N-body is due to the classical interactions with the 
other bodies of the Solar System, that labelled with GE is due to the 
post-Newtonian gravitoelectric Schwarzschild component of the solar grav- 
itational field, that labelled with J2 is due to the classical effect of the 
Sun's quadrupole mass moment J2 and that labelled with LT is due to 
the post-Newtonian gravitomagnetic Lense-Thirring component of the so- 
lar gravitational field. The LS perihelion secular rates for all the plan- 
ets, by neglecting the eccentricities of the orbits, amount to 3.95 x 10~^ " 
cy~^. For the numerical calculations of the planetary N-body precessions 
see http://ssd.jpl. nasa.gov/elem_planets. html T^ rates, For the Sun's proper 



angular momentum J, which is the source of the gravitomagnetic field, the 
value 1.9 x 10^^ kg m^ s"^ ESI has been adopted. For the Sun-planetary 
data see pi] . 



Secular rate 


Mercury 


Venus 


Mars 


'i'N-body 


1.019036 X 10^ 


8.87652 X 10^ 


2.580836 X 103 


^N-body 


-4.4630 X 10^ 


-9.9689 X 102 


-1.02019 X 103 


I^GE 


4.2981 X 10^ 


8.624 


1.351 













^J2 


5.0609 X 10-2 


5.223 X 10-3 


3.92 X 10-"^ 




-2.5375 X 10-2 


-2.613 X 10-3 


-1.96 X 10-^ 


wlt 


-3.018 X 10-3 


-4.31 X 10"^ 


-4.5 X 10-5 




1.008 X 10-3 


1.44 X 10-^ 


1.5 X 10-5 
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Table 5: First row: Formal standard deviations, in " cy~ , in the measure- 
ment of the post-Newtonian secular rates of the planetary perihelia from 
the processing of the ranging observations (E.V. Pitjeva, private communi- 
cation, 2004). They have been obtained by including also the perturbing 
effects of the asteroids which are particularly important for Mars (see Ap- 
pendixlHI). The uncertainty in the perihelion advances is one-to-three orders 
of magnitude larger than the LS effect for Mercury and Venus, respectively, 
but it is slightly smaller than the LS precession for Mars. Realistic errors 
may be an order of magnitude larger. Second row: Formal standard de- 
viations, in " cy~^, in the numerical propagation of the planetary nodal 
rates averaged over 200 years from the DE405 ephemerides (E.M. Standish, 
private communication, 2004). 



Secular rate 


Mercury 


Venus 


Mars 


u 


8.6 X 10-3 


1.037 X IQ-^ 


1 X 10^*^ 


n 


1.82 X 10"^ 


6 X 10"^ 


1 X 10-6 



masses are not known with great accuracy |29) : according to recent esti- 
mations, their mismodelled effect would amount to ~ 10"^ " cy-^for the 
secular terms (see Appendix . 

In regard to the post-Newtonian relativistic perihelion advance, it is cur- 
rently known at a 10"^ level of relative accuracy from the PPN parameters 
P and 7 [22] • It turns out that the accuracy with which the post-Newtonian 
general relativistic gravitoelectric perihelion advances can be determined is 
one-to-three orders of magnitude larger than the LS effect for Mercury and 
Venus, respectively, but it is slightly smaller than the LS precession for Mars. 
The figures of Table El for the perihelia have been obtained from the pro- 
cessing of about 300000 position observations (1911-2003) of different types: 
also the perturbing effects of the asteroids (22], which are particularly im- 
portant for Mars (see Appendix ISJ, have been accounted for. However, it 
must be noted that the realistic errors in such figures may be one order of 
magnitude larger, so that the mismodelling in the relativistic effect would 
tend to mask the LS rate. 

The mismodelled precessions induced by J2, by assuming a relative un- 
certainty of 10% in it ||22j, would be larger than the LS effect for Mercury 
and Venus; for Mars such competing classical feature would, instead, be 
smaller than the LS effect by one order of magnitude. 
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These evaluations would suggest that in a more or less near future it 
might, perhaps, be possible to analyze the perihelion of Mars only. However, 
note that a sort of a priori 'imprint' of the LS signature in the recovery of 
all the other effects affecting the perihelion of Mars would be present: it 
might make unreliable the proposed measurement by driving the outcome 
of the test just towards the desired result. Then, a better approach would 
consist in disentangling the LS effect from all the other competing effects. 
This can be done by suitably combining the perihelion of Mars with other 
Keplerian orbital elements affected by such Newtonian and post-Newtonian 
precessions. A similar approach has recently been adopted, in a Solar System 
context, in |3fl|l31j. By writing down the rates of the perihelia and the nodes 
of Mars and Mercury in terms of the LS, general relativistic, J2 and N- 
Body precessions it is possible to consider them as a linear nonhomogeneous 
algebraic system of four equations in four unknowns and to solve it with 
respect to the LS effect. The so obtained linear combination of the residuals 
of the rates of the perihelia and the nodes of Mercury and Mars is, by 
construction, independent of the effects induced by the Newtonian N-body 
and solar precessions and by the general relativistic gravitoelectric force. 
The result is 

where 

ci = 2.5275, C2 = -3.14 x 10~^ C3 = -6.67 x 10"^ (5) 

and 

Xdv = =F3.87 X 10"^ " cy-^. (6) 

It is important to note that the impact of Mercury is relatively small, as 
can be inferred from eq.©- Note that eq.Q does not cancel out the Lense- 
Thirring effect: however, according to Table |1J its nominal impact amounts 
to 2 X 10~^ " cy~^. According to the results for the nonsingular elements 
h, k, p, q of Table 4 of H2]) The formal 1-sigma error affecting eq.Q can 
be evaluated as 41%. 

5 Conclusions 

In this paper we have investigated the possibility of testing the multidi- 
mensional brane-world scenario proposed by Dvali, Gabadadze and Porrati 
by analyzing in the weak-field regime of the Solar System the pericentre 
advance derived by Lue and Starkman. 
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It has been shown that the systematic errors induced by the multipolar 
expansion of the classical part of the terrestrial gravitational potential rule 
out the possibility, suggested by Ciufolini, of using the perigee of the pro- 
posed LARES/ WEBER-SAT satellite in order to test the Dvali, Gabadadze 
and Porrati gravity by measuring the Lue-Starkman pericentre precession. 
Indeed, according to the EIGEN-CGOIC Earth gravity model from the dedi- 
cated CHAMP and GRACE missions, the errors induced by the mismodelled 
even zonal harmonics of the multipolar expansion of the Earth's gravita- 
tional potential would be up to three orders of magnitude larger than the 
LS precession of interest which amounts to ^4 x 10~'^ mas yr~^ only. More- 
over, also the observational accuracy, evaluated with optimistic assumptions, 
would not be sufficient to measure such an effect with the LARES perigee. 

On the contrary, from the ranging data to Mars and Mercury it might 
be possible, in a near future, to try to test the LS precession, which amounts 
to =f4 X 10^"*^ mas cy~^, with a suitable linear combination of the residuals 
of the perihelia and the nodes of Mars and Mercury. Indeed, the perihelion 
advance of Mars, which is the currently best tracked planet of the Solar 
System, has recently been measured with a formal standard deviation of 
1.5 X 10^^ mas over a time span of 92 years. The main systematic errors 
of gravitational origin, which severely affects a single perihelion rate, are 
instead cancelled out by the proposed linear combination. The obtainable 
accuracy could roughly be evaluated to be of the order of 41% at 1-sigma. 

6 Appendix: The impact of the asteroids on the 
Mars orbital motion 

The bodies of the asteroid's belt induce both secular and long-period pertur- 
bations on the orbital elements of Mars whose amplitudes are proportional, 
among other things, to their masses. Since they are relatively poor deter- 
mined with respect to those of the major planets of the Solar System the 
mismodelled part of the orbital perturbations induced by them might be- 
come relevant also because their nominal size is larger than the LS effect. 
It is important to note that the periods of some of the most powerful per- 
turbations amount to tens and even hundreds of years: it would be a major 
limiting factor for analysis covering relatively short time spans of available 
data. 

In the framework of the linear Lagrange-Laplace perturbation theory the 
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equation for the rate of the longitude of the perihehon is 



dw Vl - e2 &R tan(i/2) dlZ 

= 1 ^ ' ' {1\ 

dt no?e de na^Vl^ di ' ^ ^ 

where TZ is the disturbing function which accounts for the effect of a third 
body, whose parameters are denoted with on a given planet, whose pa- 
rameters are denoted without '. It is 

n = GM' (-^-^-^]. (8) 
\ |r — r I r / 

Eq.Q can be expanded as a multiple Fourier series with respect to the six 
orbital parameters of the perturbed and perturbing bodies X,zu,^l, X ,zu , 
where A is the mean longitude 

71= ^ Fjj/i^i^^/ cos{jX + j X +lw + 1 w +mft + m^}). (9) 

jj II mm 

In turn, the coefficients Fjj'u'mm'^ which depend on a, a', e, e', can 
be expanded as a Taylor series with respect to e, e', sini, sini' due to 
the smallness of the eccentricities and the inclinations of the Solar System 
bodies 



oo oo oo oo 



= E E E E^-..3..el'l--^e''"""%ml-l--3,,,|..V-.,'. 

ri=0 T2=0 r'3=0 r4=0 

(10) 

In it the coefficients Cr^r2r3r4 depend on a and a . In eq.Q the so called 
short-period term are those which contain the mean longitudes A and A , 
while the secular terms are those which do not contain them, i.e. for which 
j = j = 0. In regard to the latter ones, it turns out that {7Z), averaged 
over the mean longitudes, only contain the terms 1, e^, e , sin'^i, sin'^z', 
ee cos(ro — w ), sin i sin z cos(il — $7 ). 

Let us now specialize these general concepts to the Mars-asteroids sys- 
tem. 



According to [22]) the nominal amplitude of the secular perturbation 
on the longitude of the perihelion of Mars induced by the major aster- 
oids amount to 3 x 10~^ " cy~^. The three most massive bodies (1) 
Ceres, (2) Pallas and (4) Vesta are responsible of roughly 80% of such 
perturbation, i.e. ~ 2.4 x 10~^ " cy~^. The latest reported measure- 
ments of the masses of the asteroids (221 precise up to ~ 0.07—0.3%, 
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so that the mismodelled part of the Mars perihehon secular precession 
due to the three major asteroids amounts to ~ 1 x 10~^ " cy~^. The 
masses of the smaller asteroids like (3) Juno, (7) Iris and (324) Bam- 
berga are known to a ~ 1% level, so that the mismodelled part of their 
perturbation may be evaluated as ~ 2 x 10~^ " cy~^. The impact of 
the asteroid ring, i.e. the ensemble of the minor asteroids which can 
be modelled as due to a solid ring in the ecliptic plane [33]) can be 
worked out, e.g., with the Lagrangian approach and the disturbing 
function of the Appendix of _33 . By using the values of [221 for the 
ring's radius and mass it turns out that the secular perturbation on 
tijMars amounts to 2.8 • 10"'^ " cy~^, with an uncertainty of 3 • 10~^ " 
cy"^ 

• In regard to the short-period perturbations, according to the 
terms of eq.Q containing also the mean longitudes are, in principle, 
not negligible because of the occurrence of many orbital resonances 
which induce harmonic perturbations with periods tens or even hun- 
dreds of years long. The focus in [22] is on the perturbations of the 
mean longitudes: the nominal values of their amplitudes are of the or- 
der of 10~^ — 10"'^ " cy~^(see Table IV of |H2|). The current accuracy 
in our knowledge of the asteroids masses makes the mismodelled part 
of such perturbations smaller than the investigated braneworld effect: 
moreover, it should also be pointed out that the perturbations in the 
mean longitudes, which depend on the square of their periods, are 
larger than those in the other orbital elements, which, instead, depend 
linearly on their periods. These considerations would suggest that the 
resonant harmonic perturbations on the Martian perihelion due to the 
asteroids should not compromise the recovery of the LS effect. How- 
ever, a careful and more detailed analysis of such an important topic 
is needed. 

Concerning the error in the observational determination of the Martian 
perihelion rate, according 132] in the DE 118 and DE 200 JPL ephemerides 
it was 6 X 10"3 " cy"^ with the recent EPM2004 lAA (RAS) ephemerides 
it amounts to 1 x 10~^ " cy~^(E.V. Pitjeva, personal communication 2004). 
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